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Abstract. Results from a database analysis of H-mode electron density profiles on 
the Tokamak a Configuration Variable (TCV) under stationary conditions show that 
the logarithmic electron density gradient increases with collisionality. By contrast, 
usual observations of H-modes showed that the electron density profiles tend to flatten 
with increasing collisionality. In this work it is reinforced that the role of collisionality 
alone, depending on the parameter regime, can be rather weak and in these, dominantly 
electron heated TCV cases, the electron density gradient is tailored by the underlying 
turbulence regime, which is mostly determined by the ratio of the electron to ion 
temperature and that of their gradients. Additionally, mostly in ohmic plasmas, the 
Ware-pinch can significantly contribute to the density peaking. Qualitative agreement 
between the predicted density peaking by quasi-linear gyrokinetic simulations and 
the experimental results is found. Quantitative comparison would necessitate ion 
temperature measurements, which are lacking in the considered experimental dataset. 
However, the simulation results show that it is the combination of several effects that 
influences the density peaking in TCV H-mode plasmas. 



1. Introduction 

In recent years numerous works have been devoted to particle transport studies (see the 
most recent review [1] and references therein). The experimental observations and the 
theoretical predictions show a good overall agreement. It is now widely accepted that the 
electron density profile in the core of tokamak plasmas is mainly tailored by turbulent 
mechanisms. The turbulent state present in the plasma is intimately related to the 
density peaking and its dependence on the plasma parameters. In this work we apply 
a unified quasilinear theory of ion temperature gradient (ITG) and trapped electron 
(TEM) microinstabilities and turbulence [2] that has emerged in the last years to a very 
specific problem, namely the density profile behaviour in TCV core plasmas. The success 
of quasi-linear theory in explaining particle transport in TCV L-mode and elTB plasmas 
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[2j 13] motivates us to test the theoretical predictions against an H-mode dataset as well, 
using the same methodology. Our goal is to reproduce the observed experimental trends 
with numerical simulations and understand which transport mechanisms are influential 
on the density peaking. 

Alcator C-Mod internal transport barriers (ITBs) with very strongly peaked density 
profiles showed density profile flattening with modest electron heating, which was 
explained by the onset of TEM turbulence, using linear and nonlinear gyrokinetic 
simulations [I]. The density peaking in the C-Mod ITBs was thought to result from 
the Ware pinch, with the ITB density gradient during electron heating determined by 
the balance of Ware pinch and TEM turbulent fluxes, thus acquiring a strong inverse 
dependence on temperature. In TCV electron ITBs, where transport barriers are seen 
mainly on the electron temperature profile [3] , the density peaking results from a strong 
thermo diffusive pinch [6] component which can be explained [3] using the quasi-linear 
model that shall be used in this paper . 

The first experiments on TCV with third-harmonic X-mode electron cyclotron 
resonance heating (X3 ECH) [7J [8] reported reduction of peakedness of the density 
profile when intense electron heating was applied. A preliminary study [H] of these 
plasmas concluded that the density peaking is mainly caused by turbulent processes 
and suggested that the intense electron heating being favourable for TEM destabilization 
can possibly lead to profile flattening. In recent experimental campaigns of TCV, new 
data have been collected on ohmic (OH) and ECH H-modes, repeating the very first 
experiments and also exploring a wide parameter range. These experiments, where 
special attention was made to the quality of the collected data, allow us to better 
explore the experimental dependencies and to have more confidence in validation of the 
theory. The experimental dataset is used to test the quasi-linear model proposed in [2] . 
Our work is complementary to the recent study with AUG data [10J. 

In the next section a brief overview of the TCV H-mode experiments is given with 
a particular attention to the typical parameter ranges within which these plasmas are 
sustained. In section [3] the theoretical model that is used for the interpretation of the 
experiments is summarized. In section [4] the results of the simulations are discussed and 
compared to the experimental observations. Finally, in section [5j general discussion and 
concluding remarks are given. 

2. TCV H-mode plasmas 

In TCV (R = 0.88 m, a = 0.25 m, I p ^ lMA,B t ^ 1.5 T, k < 2.8), L-mode to H-mode 
transitions are obtained with and without additional heating [TT) IT2]. The additional 
heating is provided by powerful third-harmonic electron cyclotron heating system (X3 
ECH, 118 GHz, 3 x 0.5 MW, 2 s, top launch) [7]. In this work we focus on standard 
ELMy H-modes with type I and type III ELMs with typical ELM frequencies around 
35-100 Hz. Time traces of a typical pulse are shown in figure [TJ More 'exotic' scenarios 
such as the quiescent ELM-free [8] and snowflake [T3] H-mode scenarios are beyond the 
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Table 1. Typical parameters of the TCV H-modes. 

scope of this present paper and not considered in the analysis. 

The experimental database is built on a representative set of sufficiently diagnosed 
H-mode deuterium pulses. Some typical parameters of these plasmas in the database are 
shown in Table [T] The dataset contains pulses with OH only and those with additional 
X3 ECH, the different power levels are relatively well covered. The dataset is somewhat 
less extended in terms of changes in plasma current and contains plasmas with very 
similar shape. 

The majority ion species is deuterium and the most abundant impurity is carbon, 
originating from the wall tiles. The typical carbon concentration is a few per cent of 
the electron density. The effective charge of these plasmas, estimated from soft X-ray 
measurements [14] , is found to be Z eS ps 2 for all the pulses considered in the analysis, 
independently on the heating scheme applied. The calculated total (OH plus bootstrap) 
current, assuming neoclassical conductivity with the measured Z eS value, is consistent 
with the magnetic measurements. 

2.1. Density and electron temperature measurements 

The electron density and temperature profiles are measured via Thomson scattering. 
Density profile measurements are cross-calibrated with the far infrared interferometer 
(FIR). Typical temperature and density profiles are depicted in figure [2] The radial 
flux label throughout this paper is = yfy — ipo) / (fpb — V'o) > where ip is the poloidal 
flux, ipo and ipb are the poloidal flux at the magnetic axis and the plasma boundary, 
respectively. The derivatives of the profiles with respect to p^ are determined from cubic 
spline fits. The fit was performed on profiles collected during a stationary phase of pulses, 
i.e. where the main plasma parameters (plasma current, line averaged density, internal 
inductance) do not change significantly. The normalized inverse scale length (normalized 
gradient) of a profile is defined as: R/L x = — R (|Vp^|) dlogX/dp^, where R is the 
major radius and the (|V/ty|) is the Jacobian of the r — > p^ coordinate transformation 
with r being the distance from the magnetic axis on the outboard midplane. The profiles 
of the normalized gradients are shown in figure [3} 

All the considered pulses exhibited sawteeth. This results in flat average profiles 
within p^ ps 0.4 — 0.6 surface, depending on I p . 
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2.2. Ion temperature measurements 

The C 6+ ion temperature is measured by charge-exchange recombination spectroscopy 
(CXRS) [15]. Neutral particle analysers (NPA) [16] provide estimate of the central ion 
temperature. In these H-mode plasmas, however, due to geometrical constraints, the ion 
temperature measurements are limited outside p^ ~ 0.8. This fact, strictly speaking, 
does not allow us to measure the ion temperature profile in the core. 

figure|2]b) also shows a typical ion temperature profile between p^ = 0.8 — 1 overlaid 
on the electron temperature profiles. The ion temperature does not change significantly 
when the ECH is applied. From these typical profiles the local value of T e /T ; at p^ = 0.7 
(Pvoi ~ 0.6) is around 1 for phases with ohmic heating only, and about 2 for phases with 
ECH. 

The present ion temperature data suggests that the ion temperature gradient is 
lower than that of the electrons. Assuming the functional form T\{p) (T e /T e (p)) l/ for 
the ion temperature profile, the relation R/Lj>. = vR/Lt s follows. In figure [2] b) the 
dashed curve, derived from the electron temperature profile of #40089 using p = 0.8 
and v = 0.5, follows rather well the ion temperature measurements outside p^ ~ 0.8. 
Following these considerations we estimate R/L^ to be around 3-5 at p^ = 0.7. 

2.3. Collisionality dependence 

Statistical analysis of JET [T7] and AUG [IB] showed that collisionality is the most 
important scaling parameter for density peaking in H-modes. Both machines report 
the flattening of the density profile with increasing collisionality. Note that in these 
machines collisionality is controlled by the combination of neutral beam injection (NBI) 
heating and ECH. Despite its complexity, collisionality is a frequently used parameter to 
map parameter dependencies and for inter-machine comparisons. In this subsection the 
collisionality dependence of the density and temperature profile peaking is presented. 

Later (in section [4]), the results of the numerical modelling will be presented as a 
function of the following definition of collisionality: 

v = u ei /(v t h/R), (1) 

where u ei is the electron-ion collision frequency, v t h = ^/T-Jm the thermal velocity and 
R the major radius. In order to facilitate the comparison with previous results the 
experimental dependences are also plotted as a function of an effective collisionality 
adopted from [19]: u eS = 0.1n e Z efi R/T^(^ z/ ci /w Dc ), where Tig ; £1X6 the electron 
density in 10 19 m~ 3 and the electron temperature in keV with Z c g = 2. 

figure [4] shows the normalized logarithmic gradients of the electron temperature 
(R/LtJ and density (R/L n ) profiles at p^ = 0.7 as a function of collisionality v and v^. 
Points around v ~ 0.02, the most typical TCV OH H-mode pulses, have R/L n around 
3. Samples from reversed field campaigns with unfavourable ion V-B drift direction with 
somewhat higher density, are located at around v = 0.03 having R/L n rs 3.5. In the 
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0.008 < v < 0.02 range, which is reached by adding 0.5-1.5 MW ECH, R/L n decreases 
to about 1.5. 

Note that the global trend of R/L n as a function of the collisionality is very similar 
to the one observed on FTU fully non-inductive, high density electron heated plasmas 
[2D] . By contrast, H-modes at JET [T7] and AUG [TD] showed that R/L n is decreasing 
with increasing collisionality. 

3. Theoretical model 

In Ref. [2] a quasi-linear gyrokinetic framework has been proposed for particle transport 
with which TCV L-mode plasmas and electron internal transport barriers [3] were 
studied extensively. In this paper we strive to test this model against TCV H-mode data 
to understand the density profile behaviour. Although there is no general agreement 
in the community on the applicability of different quasi-linear models, significant effort 
has been made to compare various models against non-linear simulations [2TI [22l |23[ 121] 
and good agreement was found. It was also shown that the fully developed turbulence 
preserves many features of the linear evolution [22], therefore with an ad-hoc weighting 
of the linear spectra, non-linear fluxes can be acceptably reproduced [23]. Moreover, 
quasi-linear predictions fit very well to experimental observations [21 |3j [TD]. A relevant 
exception to this occurs near marginal stability, where a nonlinear upshift of the TEM 
critical density gradient has been found [I], which increases with collisionality [25]. This 
effect is clearly outside the framework of quasilinear theory. EC heated experiments are 
in general not near marginal stability to temperature gradient driven TEMs [2H] with 
R/Lt,, ~ 8 [2Z] but may be near marginal stability to longer wavelength density gradient 
driven TEMs [I]. In our dataset, R/L n ranges from about 1 to approximately 3.5, with 
the upper value lying above the density gradient driven TEM treshold (as shown later 
in figure [5]) . 

We solve the linear gyrokinetic (GK) equation for electrostatic perturbations with 
a wavenumber k = (0, k y ) to obtain w*, 1 k , the real and imaginary part of the 
mode frequency and the particle flux, respectively, of the most unstable mode, cj k > 
represents a mode turning in the ion diamagnetic (ITG) direction, while u k < 
corresponds to the electron diamagnetic direction (TEM). 

, 7 fc , T fc are evaluated on a range of k y . We then obtain the quasi-linear fluxes 
and the average mode frequency of the turbulent state by a weighted average over the 
mode spectrum: 



where R can stand for T k , q k , u k . The Wk weights are usually chosen according to a 
quasi-linear rule of the form: 




(2) 




(3) 
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We use the same parameters £ = 2, Aq — 1 as in |2]. We shall discuss effect of this 



choice later, in section 4.6 



In order to separate the different mechanisms that are responsible for the density 
peaking, one can decompose the particle flux such as: 

T = A^ + B-^ + C, (4) 

Lin -^Te 

which, assuming stationary plasma and using the zero flux condition when no particle 
sources are present (as is the case in the core of these ECH/OH plasmas), transforms 
into 

R 

Ln 



R 

-Ctj Cp- (5) 



The thermodiffusion coefficient Cr = —B/A and the other pinch coefficient Cp = —C/A 
are evaluated numerically ^ 12] • 

4. Simulation results 

4-1. Setup for the simulations 

Linear gyrokinetic calculations have been performed with the initial value flux-tube 
code GS2 [28]. We introduce two model cases representing TCV OH and ECH H-mode 
pulses. For both R/Lt b = 9, R/L^ = 6, Z c s = 2. For the OH reference case T e /T; = 1.0, 



v = 0.024, for the ECH case TJT-, = 2.0, v = 0.008 (c. f. subsection |2TT[) . 

The magnetic equilibrium of a representative discharge (see figure [T| was calculated 
by the equilibrium code CHEASE [29|. The = 0.7 flux surface with q = 1.2, s = 0.7 
was used for flux-tube simulations. For the collisions both pitch angle scattering and 
energy diffusion are taken into account. 13 k y values are used in the range of [0.08, 1.6] 
distributed in a logarithmic way. All growth rates and frequencies are normalized to 
v[ h /R. The binormal wavenumber k y p\ is normalized to pi, the normal wave number 
is always zero (k x = 0). Note that on the outboard midplane the normal direction is 
radial. The timestep was fixed to At = 0.05 in units of R/v\ h . 32 grid points for each 
2tt turn in 6, and 12 poloidal periods were used. The parameters are then scanned 
according to the various simulation results presented in the next subsections. 

4-2. Finding the stationary density gradient 

In the absence of core particle sources, as it is the case in the TCV OH/ECH H-modes 
[30], the stationary value of the density gradient R/L n \ sUt corresponds to the T = 
condition. Here we only consider the particle flux ascribed to turbulence, we shall deal 



with neoclassical contributions later, in section 4.5. Other contributions to the total 
particle flux (e.g. MHD, ripple losses, etc) are neglected. It is also assumed that these 
approximations hold under both OH and ECH conditions. R/L n \ tat is obtained by 
making a set of simulations for different values for R/L n , keeping all the other input 
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parameters fixed. At each R/L n , (T) is evaluated with equation [2] Then the zero point 
is found by linear interpolation. 

In order to understand the relationship between the linear and quasi-linear fluxes, 
in figure [5] we show the different contributions in equation [2] together with the mode 
frequency and growth rate of the most unstable mode as a function k y pi for the OH 
reference case. First we note that the quasi-linear rule gives small weight at large 
wavenumbers (short scales), therefore (r) and (u r ) are mainly determined by the 
kyPi < 0.4-0.5 wavenumber range (ITG and TEM range). 

When the density gradient R/L n is small (figure [5] a), b)), the most unstable 
modes are ITGs (positive u r ), the WkT k contributions to (r) throughout the covered 
wavenumber range is negative (inwards), making (r) pointing unambiguously inwards. 
Increasing R/L n (figure [5] c) and d)), the growth rate of the most unstable modes 
decreases and at low k y pi, WkT k > contributions appear, making (r) = when 
R/L n = R/L n \ Btat is reached. Further increasing R/L n ((figure [5] e) and f)) the mode 
frequency changes sign and TEM modes appear at small wavenumbers. When these 
modes dominate, the particle flux (r) is positive (outwards). 

figure [6] shows the quasi-linear particle flux (T) normalized to the total electron 
heat flux as a function of R/L n for the two reference cases. Each point in the figure is 
a result of 13 simulations (averaged over the k y spectrum as in figure [5]). The average 
mode frequency (u T ) is also plotted showing that, indeed, the zero flux point is found 
between ITG and TEM type modes [2] . Note the role of ITG-TEM balance at a certain 
density gradient in the spectrum and also in adjusting the value of R/L n \ stat . It must 
also be stressed out, that at the co-existence of different modes, one must certainly 
account for non-linear effects jH I25J [2U EE]- It is not expected that the absolute value 
of i?/L n | stat from the quasi-linear estimate of the particle flux instantly matches the 
experimental values; however, its parametric dependence on the plasma parameters is 
remarkably well retained (2J El [TO] - 



4-3. Veg-T e /Ti scan 



We endeavour to reproduce the collisionality dependence reported in Sec. |2.3| with 
numerical simulations. Since in TCV H-mode experiments collisionality is mainly 
controlled by intense electron heating, the T e /Ti ratio certainly increases with additional 
ECH. It is expected that the temperature ratio has a significant effect on the nature of 
the dominant instabilities and hence on the particle flux. We also want to explore the 
effect of the electron-ion collisions at each considered temperature ratio. Therefore, we 
perform a double parameter scan in the range of v = [0.008, 0.012, 0.016, 0.020, 0.024] 
together with a scan in T e /1] = [1, 1.5, 2]. 

The results are summarized in figure [7] a), which shows the predicted value of the 
density gradient R/L n \ stat as a function of the collisionality v. Each point is obtained 
from a scan over R/L n to determine R/L n \ stat where (r) = 0. Therefore each point in 
figure[7]is the result of about 10 x 13 linear simulations. When T e /T; = 1.0 the predicted 
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density gradient is decreasing with collisionality. As the temperature ratio is increased, 
this trend is reversed and at T e /Tj = 2.0 an increase of R/L n \ stat is predicted with 
increasing v. These results show a clearer trend (figure [7] b)) if now one plots R/L n \ stat 
as a function of (u T ) which characterizes the nature of the background turbulence and 
orders very well the predicted values of _R/L n | stat [21 HD] and also clearly shows that the 
change in the dominant instabilities alters the collisionality dependence. Therefore, an 
explanation of the different behaviour of TCV H-modes with respect to JET/ AUG 
plasmas can be due to the different T e /T[ ratio and hence a different collisionality 
dependence. 

Collisionality provides a purely convective term to the particle flux pQ which is 
directed outwards for ITG and inwards for TEM modes and increases with increasing 
collisionality. Intuitively one expects that increasing v tends to push u T towards more 
positive values p] (towards ITG type turbulence) due the stabilizing effect of collisions 
on TEM turbulence [H [32] . The net effect of collisionality on R/L n \ t . stems from 
these two effects: v influences the turbulence regime and, depending on the turbulence 
regime, the particle flux. Since an increase in collisionality implies an increase in the 
real mode frequency the effect of the collisionality on the TEM modes driven outward 
flux is expected to be weak pQ. Indeed, in figure [7] we observe almost no collisionality 
dependence at T e /Ti = 1.5, where (u r ) ~ 0. Moving away from this point toward 
either more positive or more negative (u T ), the net particle flux is less inwards, and 
R/L n \ stat is found at a lower gradient. The increase of T e /T; at fixed T ; is destabilizing 
for ITG and enhances the TEM [33] activity as well. Their interplay is such that TEM 
takes over rapidly with increasing T e /T\, in particular with ECH i.e. by increasing T e 
[34]. Indeed, starting from the OH reference point in the ITG instability domain, with 
increasing T e /Ti and decreasing z>, (u r ) decreases and the predicted value of R/L n \ stat 
is increasing. In the transition region between ITG and TEM, the largest R/L n \ stat is 
obtained, then moving towards the TEM domain the density peaking decreases, leaving 
the point corresponding to the ECH reference case with a very similar R/L n \ stat value. 

We conclude that albeit ECH has a significant effect on the values of v and T e /T- U 
the change in these parameters only is not sufficient to explain the observed overall 
experimental behaviour. 



4-4- R/Lt b andR/L Ti dependence 

It is expected that the temperature gradients have a large effect on the predicted value 



of the density gradient. As it was pointed out in Sec. |2.2[ the ion temperature gradient 
is generally lower than that of the electrons even at the high collisionality OH plasmas. 
The ions are even more decoupled from the electrons when intense electron heating is 
applied. Thus we expect that R/L^ does not change significantly in these plasmas 
(more precisely it does not increase). The experimental data are more clear on the 
R/L Te dependence; however it is useful to study the dependencies on the parameter, as 
well. 
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Two scans of Rj Lt £ in the range of [6, 7.5, 9, 10.5, 12] have been performed: the first 
with OH only parameter set (T e /T; = 1.0, v = 0.024); the second with ECH parameter 
set (T e /Tj = 2.0, v = 0.008). The other parameters were kept fixed (c. f. Sec. 4.1). In 
the same way, we change the value of R/L^ in the range of [4, 5, 6, 7.5, 9]. 

In figure [8] we present the results of these four scans in terms of 7] ei = 
(R/Lt £ ) I (R/Lxi)- Figure [8] a) shows that the change in R/L^ has a different effect 
on the OH and the ECH reference cases. Indeed, starting with the OH reference case 
(triangles), increasing R/L^ predicts smaller R/L n . On the other hand, this parameter 
has the opposite effect on the predicted peaking for the ECH case. In figure [8] b), it is 
shown that increasing R/Lx e increases R/L n for the OH case, while the density peaking 
is rather insensitive to the change in the electron temperature gradient for the ECH 
case. 

In figure [9] we decompose the predicted R/L n in terms of Cr and Cp as in equation 
(|5| in order to identify the change in the thermodiffusive and the other contributions. We 
see that the larger contribution is always the thermodiffusive pinch CtR/Lt,, and usually 
the two pinches change in the opposite direction and a complicated interplay between 
the two yields the predicted R/L n . This means that the parametric dependence of the 
density peaking is largely determined by the starting reference parameters as mentioned 
above. 



In figure 10 a) we plot these results against (ui T ) which, again, orders very well 
the simulation results R/L n [21 EH]- As it was already shown in Sec. 4.3, the OH 
reference (shaded triangle) case is in the ITG regime ((u; r ) > 0), while the ECH 
(shaded square) reference case is predicted to be more TEM ((u t ) < 0). Increasing 
the electron temperature gradient results in moving to the left on the (u T ) axis (towards 
negative values), while the ion temperature gradient pushes towards larger (u T ) values. 
The density peaking is maximal around (u T ) ~ and the predicted value decreases 
towards both ITG and TEM regime as predicted for elTB and L-mode parameters 
[31 [2]. We therefore show that this is a very general feature. In figure 10 b) c) we see 
the thermodiffusive contribution increases towards the ITG regime, while the Cp term 
decreases. This behaviour of the two pinch terms reinforces their suggested universal 
role in particle transport [2]. 

As mentioned in section 2.2 the ion temperature measurements indicate that 
(R/L T .)/ (RJ L Te ) « 0.5, although we do not have measurements near p^ = 0.7, where 
we have based our analysis. If this ratio is further decreased in the simulations, 



moving left in figure 10 a), the value of i?/L n | stat decreases, following the trend seen 
previously. However the range of h y p\ where the ITG is the most unstable mode, 
typically k y pi = 0.1 — 0.6, shrinks and slab-like electronic modes very elongated along 
the magnetic field lines appear as seen in Refs. [351 Although their contribution 
to the particle flux is small, so is that of the remaining ITG and TEM. This domain 
requires non-linear simulations and the role of non-adiabatic passing electrons to be 
studied in detail, which is left for further analyses. 
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4-5. Effect of the Ware-pinch 

In ohmic plasmas, the neoclassical Ware-pinch can have significant contribution to the 
density peaking jH [37] . The strength of this effect depends also on the heating scheme 
applied [U HO]- In the presence of the Ware-pinch, one needs to find R/L n \ sUt from: 
T + n e W p = [1]. In order to cancel normalization factors, it is more convenient to use 
the T/Q e ratio; therefore, we divide with the electron heat flux Q e : 

T ^ n c W p q . . 

Qe Qe 

We evaluate the first term using the ratio of the normalized fluxes (r) / (Q e ) (which is 
proportional to T e T/Q e ) from numerical simulations, while the second term is computed 
from experimental measurements using Q e = Qg Xp and power balance (PB) arguments 
assuming a stationary plasma, where 

Qr P = -Xe PB -eVT e = Xe PB ^^- (7) 



We find 



(r> rw p 



from which one can find the stationary density gradient value with the Ware-pinch effect 
included. Note that W p < 0. Note also that for any consistent normalization of (r) and 
(Q e ) , T e of equation [7] cancels in equation [8j 

We compare the Ware-pinch contribution for the OH and ECH reference cases 
substituting typical values of W p and Xe B i n equation [8j assuming that the intense ECH 
does not directly drive significant particle flux that must be countered by the turbulent 



flux to restore the total flux to zero. In figure 11 the ratio of the normalized fluxes 
(r) / (Q c ) without (symbols), and with (no symbols) the Ware-pinch contribution are 
shown as a function of R/L n ; (a) shows the OH reference case with = 0.3m 2 s -1 , (b) 
is the ECH reference case with Xe B = lm 2 s -1 . For both cases W p = — 0.3ms _1 . It can 
be seen that the shape of the (r) curve around R/L n \ stat can influence the position of the 
new value of R/L n \ staV The addition of the Ware-pinch results in a vertical downshift of 
the curve. At the OH reference case (r) is almost linear with R/L n . Due to its modest 
slope, the upshift in R/L n \ stat can be as large as 1. For the ECH case the slope of the T 
curve is somewhat reduced which would imply a larger change in i?/L n | stat . However, 
for the ECH case the contribution of the Ware-pinch is negligible mainly because of 
two reasons. First, the loop voltage decreases as the electron temperature increases at 
constant current [I]. This reduces the absolute value of W p . Second, due to the profile 
stiffness the value of Xe B R/-^T e is larger in ECH plasmas, which makes the contribution 
of the second term in equation [8] much smaller. 

4-6. Choice of the quasi-linear rule 

Our analysis is consistent with the observed experimental data. However it is based 
on a specific quasi-linear rule (equation [3] with £ = 2). Other similar studies reported 
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successful interpretations with different quasi-linear theory pm E31 EH [38], [22], [3] using 
different summation rule, or even using only one k y mode. A question arises whether the 
results reported here are strongly dependent on the choice of the w weights in equation 

El 

We evaluated the quasi-linear fluxes and R/L n \ st , with several different choices of 
w for the parameter scans discussed previously. In figure 12 we show again the results 
of the temperature gradient scans (see figure 10 ) summing many modes with the power 
law used thorough this paper and Ref. [2] (squares), summing with the classical mixing- 
length estimate (diamonds), summing with the exponential rule (triangles) motivated 
by fluctuation measurements [23], using one mode where the mixing length transport 
is maximal (circles) [3] and using one mode with k y pi = 0.3 (stars) [10J. Although, the 
differences in the predicted absolute value of R/L n \ stat can be as large as 1.5, the relative 
changes in i?/L n | stat and the parametric dependencies are nearly identical, independent 
on the choice of the quasi-linear rule. Figure 12 and 10 show that the main trends do 
not depend so much on the quasi-linear rule but rather on the main dependence of the 
underlying most unstable modes. However, it is safer to use the whole spectrum in 
order to be sure to capture both the effects of ITG and TEM modes. The importance 
of the most unstable modes explains why similar trends can be expected in non-linear 
simulations [221 12H EI] and in the experiments. 



5. Conclusions 



The collisionality dependence of the local logarithmic electron density gradient R/L n 
has been investigated in ohmic and electron cyclotron resonance heated H-mode plasmas 
on TCV. The density profile flattens when ECH power is added, the normalized density 
gradient R/L n \ stsct at p,^ ~ 0.7 decreases from about 3 to approximately 1.5. The 
dependence of R/L n \ stat on collisionality is rather unusual: the density peaking increases 
with increasing collisionality. With local quasi-linear gyrokinetic simulations it was 
shown that the change in collisionality alone cannot reproduce the observed experimental 
behaviour. The simulation results show that the most important parameter is the 
ion temperature gradient R/L^, which was not available in these experiments. Using 
R/Lj-. = 6 puts the OH and the ECH reference case in the ITG and the TEM regime, 
respectively, with predicted R/L n 2.5 and 3. For the OH case the contribution of 
the Ware-pinch can be as large as AR/L n rj 1 resulting in a peaking around 3.5 which 
is in good agreement with the experiments. However, the flat ECH profiles are not 
recovered. If a lower ion temperature gradient R/L^ = 4 is assumed, the predicted 
OH density gradient increases to around 4 (with the Ware pinch included) and that 
of the ECH case decreases to 2. The experimental trends are, therefore, qualitatively 
explained. The edge ion temperature data suggest that R/L T . is about 3-5, however 
new experimental data are required to compare with our simulation results. We have 
also seen that at very low R/L^, below 3-4, electron modes with large k y p\ are seen 
to play a role, therefore new simulations are also required in this domain with non- 
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adiabatic passing electrons and non-linear model. The quasi-linear simulations predict 
that R/L n \ stat for the OH and the ECH reference case decreases as Rj Lx i is decreased, 
however dedicated experiments are needed to explore how the T\ profile changes between 
OH and ECH phases. 

The simulation results, where the parameters can be changed individually, also 
show that the collisionality dependence can actually be due to a change in T e /T; and 
(R/LTe/R/L/Tj) rather than the collisionality itself. Table [2] summarizes the change in 
i?/L n | stat due to the considered parameters. 

The (r) = stationary condition is set by the balance between ITG and TEM 
modes as obtained in L-mode simulations [2J. Depending on the turbulence regime, 
the dependence of i?/L n | stat on the collisionality, temperature ratios and temperature 
gradients can be very different, but in all cases the average mode frequency (u t ) as 
a figure of merit of the background turbulence clarifies the simulation results. This 
property of real mode frequency of the most unstable modes was observed in the study 
of L- and H-mode AUG plasmas [19], \10\ , L-mode and elTB TCV plasmas [31 [2j , and now 
has been confirmed for TCV H-mode parameters, as well. The neoclassical Ware-pinch 
can also contribute to the peaking in certain conditions typically when the stationary 
point lies in an ITG dominated regime and when the electron heat flux is not too large. 

It is interesting to relate the present results to those of AUG reported recently in 
[TO] . The density peaking in those considered H-mode plasmas increases with increasing 
additional ECH power, mainly due to the increase in the ratio of the logarithmic electron 
temperature gradient to the logarithmic ion temperature gradient and partly also from 
the increase in the electron to ion temperature ratio. Figure 7. (a) in [10J can be 



directly compared to the curve with stars in figure 12. In purely NBI heated AUG 



H-mode plasmas ITG modes are the most unstable modes, positioned on the right part 
of the plot with positive frequencies. Adding ECH destabilizes the TEMs and R/L n \ t t 
increases as the mode frequency decreases. It is predicted that even larger amount 
of ECH would result in decreasing peaking [10], however, this could not be validated 
by those experimental results. According to our interpretation, the TCV OH plasmas 
are positioned very close to the ITG-TEM boundary and the additional ECH moves 



R/L n \ stat down on the left branch of the curve in figure 12 Therefore our results 
combined with AUG observations [10J recover the full non-monotonic (u T ) dependence. 

In these TCV H-mode plasmas, the local behaviour of the density profile can be well 
explained in a unified quasi-linear gyrokinetic theory of ITG and TEM instabilities and 
turbulence. This work confirms that the model used to understand particle transport 
in L-mode [2] and elTB [3] plasmas is also able to reproduce the main trends observed 
in H-modes. The results of a quasi-linear theory used in this work can be further 
improved and validated by non-linear simulations and also global codes are expected to 
predict more precise values of _R/L n | stat . For a better test of this theory, systematic and 
precise ion temperature profile measurements are desirable. The systematic study of 
the dependence of local density peaking on the plasma parameters such as collisionality, 
electron to ion temperature ratio, electron and ion temperature gradient, Ware-pinch 
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Table 2. Change in R/L n when the related parameter increases. 

show that the value of the plasma parameters from which these scans are performed 
are important in order to understand their effects on the stationary value of the density 
gradient. This complex interdependence is found to be greatly simplified when analysed 
in terms of the quasi-linear real mode frequency of the most unstable modes. 
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Figure 1. Time traces of (a) D Q emission; (b) central and peripheral FIR chords; (c) 
central electron and ion (C 6+ ) temperature; (d) OH and ECH power. The stationary 
profiles are taken from the shaded intervals 1.2-1.4 s and 1.46-1.6 s. 
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Figure 2. Profiles comparing conditions with ohmic heating only (diamonds, 
#38012 1.46 - 1.6 s ) to those with additional 350 kW (squares, #40089, 0.7 - 1.61 s), 
675 kW (triangles pointing down, #38012 1.2 - 1.4 s) ECH. (a) electron density, (b) 
electron temperature, ion temperature (points, #40089, 0.7 — 1.61 s). Solid lines show 
the piecewise polynomial fits which were used to compute the gradients. The shaded 
areas mark the confidence intervals corresponding to 95% confidence level. The dashed 
line represents 11(0.8) (T c /T c (0.8)) - 5 (#40089, 0.7- 1.61 s). 




Figure 3. Profiles of the normalized logarithmic gradients calculated from the 
piecewise polynomial fits depicted in figure [2] The shaded areas mark the confidence 
intervals corresponding to 95% confidence level. 
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Figure 4. Scatter plots of the normalized logarithmic gradients at — 0.7 as a 
function of collisionality over the considered database, (a) R/L n (b) R/Lt s - Different 
symbols indicate different heating phases, with ohmic heating only (triangles pointing 
down) and with central X3 ECH (squares). The colorscale (online) represents T o (0), 
the central value of the electron temperature. 



Understanding the core density profile in TCV H-mode plasmas 



19 




i.O 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 

kypi 




i.O 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 



1.0 
0.5 
0.0 
-0.5 
-1.0 
-1.5 
-2 




rt/<Q e > 

W It/ (Q e > 

tfdkyw r k e /{Q e ) 




l.O 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 




1.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 

kyPi 



i.O 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 

kypi 



1.5 
1.0 
0.5 
0.0 
-0.5 
-1.0 
-1.5 
-2.0 




7 



0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 

kyPi 



Figure 5. Left: particle flux T (diamonds); right: mode frequency u> and growth 
rate 7 (squares and circles, resp.) as a function of kyp^. (a), (b) low R/L n where ITG 
is dominant and T negative (inwards); (c),(d) R/L„ w R/L n \ st&t where r«0; (e), (f) 
high R/L n with TEM dominant and positive (outwards) particle flux. The simulation 



parameters correspond to the OH reference case (see section 4.1 ) 
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Figure 7. Predicted values of R/L n for different values of the temperature ratio 
T e /Tj, as a function of (a) the collisionality v and (6) the average mode frequency 
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Figure 8. Predicted values of normalized density gradient R/L n as a function of 
the temperature gradient ratio LtJLt c changing i£/£r« (solid) and R/Lt c (dashed). 
Different symbols indicate different heating phases, with OH only (triangles) and ECH 
(squares). The larger shaded symbols show the reference cases. 
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Figure 9. Thermodiffusive —CtR/Lt c and the other pinch — Cp contributions to 
the predicted value of R/L n (c. f. figure [8| as a function of the temperature gradient 
ratio LtJLt c changing R/L^ (solid) and R/Lt c (dashed). Different symbols indicate 
different heating phases, with OH only (triangles) and ECH (squares). The larger 
shaded symbols show the reference cases. 
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Figure 10. The predicted value of R/L n) the thermodiffusivc —C^R/Lt c and 
the other pinch — Cp contributions as a function of the average mode frequency (ui T ) 
changing R/L^ (solid) and R/Lt c (dashed). Different symbols indicate different 
heating phases, with OH only (triangles) and ECH (squares). The larger shaded 
symbols show the reference cases. 
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Figure 11. Normalized particle flux as a function of R/L n , with (no symbols) and 
without (diamonds) the Ware-pinch contribution according to equation ffify: (a) OH 
model scenario, = 0.3 m 2 /s (b) ECH model scenario, x^ B = 1 m 2 /s. For both 
cases: W p = —0.3 m/s. 




Figure 12. Same as figure 10 (a), using different quasi-linear rules: (a) R/Lxi scan, 
(6) R/Ltc scan. Different symbols show different prescriptions for the w weights in 
equation [2j The larger shaded symbols mark the reference cases as in figure 10 the 
OH case for (to x ) > and ECH case otherwise. 



